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The Casimir effect: some aspects 
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We start this paper with a historical survey of the Casimir effect, showing that its origin is related to ex- 
periments on colloidal chemistry. We present two methods of computing Casimir forces, namely; the global 
method introduced by Casimir, based on the idea of zero-point energy of the quantum electromagnetic field, and 
a local one, which requires the computation of the energy-momentum stress tensor of the corresponding field. 
As explicit examples, we calculate the (standard) Casimir forces between two parallel and perfectly conducting 
plates and discuss the more involved problem of a scalar field submitted to Robin boundary conditions at two 
parallel plates. A few comments are made about recent experiments that undoubtedly confirm the existence of 
this effect. Finally, we briefly discuss a few topics which are either elaborations of the Casimir effect or topics 
that are related in some way to this effect as, for example, the influence of a magnetic field on the Casimir effect 
of charged fields, magnetic properties of a confined vacuum and radiation reaction forces on non-relativistic 
moving boundaries. 



I. INTRODUCTION 



A. Some history 



The standard Casimir effect was proposed theoretically 
by the dutch physicist and humanist Hendrik Brugt Gerhard 
Casimir (1909-2000) in 1948 and consists, basically, in the 
attraction of two parallel and perfectly conducting plates lo- 
cated in vacuum UJ]. As we shall see, this effect has its origin 
in colloidal chemistry and is directly related to the dispersive 
van der Waals interaction in the retarded regime. 

The correct explanation for the non-retarded dispersive 
van der Walls interaction between two neutral but polariz- 
able atoms was possible only after quantum mechanics was 
properly established. Using a perturbative approach, London 
showed in 1930 ||2j| for the first time that the above mentioned 
interaction is given by VLo„(r) w — (3/4)(?i(aoa^)/r^, where 
a is the static polarizability of the atom, (Bq is the dominant 
transition frequency and r is the distance between the atoms. 
In the 40's, various experiments with the purpose of studying 
equilibrium in colloidal suspensions were made by Verwey 
and Overbeek isj]. Basically, two types of force used to be 
invoked to explain this equilibrium, namely: a repulsive elec- 
trostatic force between layers of charged particles adsorbed by 
the colloidal particles and the attractive London-van der Waals 
forces. 

However, the experiments performed by these authors 
showed that London's interaction was not correct for large dis- 
tances. Agreement between experimental data and theory was 
possible only if they assumed that the van der Waals interac- 
tion fell with the distance between two atoms more rapidly 
than 1/r^. They even conjectured that the reason for such a 
different behaviour for large distances was due to the retar- 
dation effects of the electromagnetic interaction (the informa- 
tion of any change or fluctuation occurred in one atom should 
spend a finite time to reach the other one). Retardation ef- 
fects must be taken into account whenever the time interval 
spent by a light signal to travel from one atom to the other 
is of the order (or greater) than atomic characteristic times 



{r/c > l/cOMnn, where co,„„ are atomic transition frequencies). 
Though this conjecture seemed to be very plausible, a rigorous 
demonstration was in order. Further, the precise expression 
of the van der Waals interaction for large distances (retarded 
regime) should be obtained. 

Motivated by the disagreement between experiments and 
theory described above, Casimir and Polder |01 considered 
for the first time, in 1948, the influence of retardation ef- 
fects on the van der Waals forces between two atoms as well 
as on the force between an atom and a perfectly conducting 
wall. These authors obtained their results after lengthy calcu- 
lations in the context of perturbative quantum electrodynam- 
ics (QED). Since Casimir and Polder's paper, retarded forces 
between atoms or molecules and walls of any kind are usu- 
ally called Casimir-Polder forces. They showed that in the 
retarded regime the van der Waals interaction potential be- 
tween two atoms is given by VRet{r) = — 23?M7aAOCB/(47tr^). 
In contrast to London's result, it falls as l/f . The change 
in the power law of the dispersive van der Waals force when 
we go from the non-retarded regime to the retarded one 
{Fnr ^ IjP Fr ^ 1/r^) was measured in an experiment 
with sheets of mica by D. Tabor and R.H.S. Winterton ^ 
only 20 years after Casimir and Polder's paper. A change 
was observed around 150A, which is the order of magnitude 
of the wavelength of the dominant transition (they worked in 
the range 50A — 300A, with an accuracy of ±4A). They also 
showed that the retarded van der Waals interaction potential 
between an atom and a perfectly conducting wall falls as 1 / r^, 
in contrast to the result obtained in the short distance regime 
(non-retarded regime), which is proportional to 1/r^ (as can 
be seen by the image method). Casimir and Polder were very 
impressed with the fact that after such a lengthy and involved 
QED calculation, the final results were extremely simple. This 
is very clear in a conversation with Niels Bohr. In Casimir' s 
own words 

In the summer or autumn 1947 (but I am not ab- 
solutely certain that it was not somewhat earlier 
or later) I mentioned my results to Niels Bohr, 
during a walk. "That is nice", he said, "That is 
something new. " I told him that I was puzzled 
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by the extremely simple form of the expressions 
for the interaction at very large distance and 
he mumbled something about zero-point energy. 
That was all, but it put me on a new track. 

Following Bohr's suggestion, Casimir re-derived the results 
obtained with Polder in a much simpler way, by computing the 
shift in the electromagnetic zero-point energy caused by the 
presence of the atoms and the walls. He presented his result 
in the Collogue sur la theorie de la liaison chimique, that took 
place at Paris in April of 1948: 

I found that calculating changes of zero-point en- 
ergy really leads to the same results as the calcu- 
lations of Polder and myself... 

A short paper containing this beautiful result was published 
in a French journal only one year later |0]. Casimir, then, 
decided to test his method, based on the variation of zero-point 
energy of the electromagnetic field caused by the interacting 
bodies in other examples. He knew that the existence of zero- 
point energy of an atomic system (a hot stuff during the years 
that followed its introduction by Planck |7|]) could be inferred 
by comparing energy levels of isotopes. But how to produce 
isotopes of the quantum vacuum? Again, in Casimir's own 
words we have the answer [8J : 

if there were two isotopes of empty space you 
could really easy confirm the existence of the 
zero-point energy. Unfortunately, or perhaps for- 
tunately, there is only one copy of empty space 
and if you cannot change the atomic distance 
then you might change the shape and that was 
the idea of the attracting plates. 

A month after the Collogue held at Paris, Casimir presented 
his seminal paper iQl] on the attraction between two parallel 
conducting plates which gave rise to the famous effect that 
since then bears his name: 

On 29 May, 1948, T presented my paper on 
the attraction between two perfectly conducting 
plates at a meeting of the Royal Netherlands 
Academy of Arts and Sciences. It was published 
in the course of the year.. 

As we shall see explicitly in the next section, Casimir ob- 
tained an attractive force between the plates whose modulus 
per unit area is given by 

M^0013 ^ ^ (1) 
L-^ {a/^im) cm^ 

where a is the separation between the plates, the area of 
each plate (presumably very large, i.e., L a). 

A direct consequence of dispersive van der Waals forces 
between two atoms or molecules is that two neutral but polar- 
izable macroscopic bodies may also interact with each other 
However, due to the so called non-additivity of van der Waals 
forces, the total interaction potential between the two bodies 



is not simply given by a pairwise integration, except for the 
case where the bodies are made of a very rarefied medium. 
In principle, the Casimir method provides a way of obtaining 
this kind of interaction potential in the retarded regime (large 
distances) without the necessity of dealing explicitly with the 
non-additivity problem. Retarded van der Waals forces are 
usually called Casimir forces. A simple example may be in 
order Consider two semi-infinite slabs made of polarizable 
material separated by a distance a, as shown in Figure 1 . 



a 




FIG. 1: Forces between molecules of the left slab and molecules of 
the right slab. 

Suppose the force exerted by a molecule A of the left slab 
on a molecule B of the right slab is given by 

fAB = — — Fab 7 
''ab 

where C and y are positive constants, r^B the distance between 
the molecules and t^b the unit vector pointing from A to B. 
Hence, by a direct integration it is straightforward to show 
that, for the case of dilute media, the force per unit area be- 
tween the slabs is attractive and with modulus given by 

Fslahs C 

Area ~ ' 

where C' is a positive constant. Observe that for y = 8, which 
corresponds to the Casimir and Polder force, we obtain a force 
between the slabs per unit area which is proportional to 1 /a^. 
Had we used the Casimir method based on zero-point energy 
to compute this force we would have obtained precisely this 
kind of dependence. Of course, the numerical coefficients 
would be different, since here we made a pairwise integration, 
neglecting the non-additivity problem. A detailed discussion 
on the identification of the Casimir energy with the sum of 
van der Waals interaction for a dilute dielectric sphere can be 
found in Milton's book |@| (see also references therein). 

In 1956, Lifshitz and collaborators developed a general the- 
ory of van der Waals forces ifioll . They derived a powerful 
expression for the force at zero temperature as well as at fi- 
nite temperature between two semi-infinite dispersive media 
characterized by well defined dielectric constants and sepa- 
rated by a slab of any other dispersive medium. They were 
able to derive and predict several results, like the variation of 
the thickness of thin superfluid helium films in a remarkable 
agreement with the experiments ifTlll . The Casimir result for 
metallic plates can be reobtained from Lifshitz formula in the 
appropriate limit. The Casimir and Polder force can also be 
inferred from this formula @] if we consider one of the media 
sufficiently dilute such that the force between the slabs may 
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be obtained by direct integration of a single atom-wall inter- 
action 111211 . 

The first experimental attempt to verify the existence of 
the Casimir effect for two parallel metallic plates was made 
by Sparnaay ifisll only ten years after Casimir's theoretical 
prediction. However, due to a very poor accuracy achieved 
in this experiment, only compatibility between experimental 
data and theory was established. One of the great difficulties 
was to maintain a perfect parallelism between the plates. Four 
decades have passed, approximately, until new experiments 
were made directly with metals. In 1997, using a torsion 
pendulum Lamoreaux lfl4ll inaugurated the new era of experi- 
ments concerning the Casimir effect. Avoiding the parallelism 
problem, he measured the Casimir force between a plate and a 
spherical lens within the proximity force approximation ifTsll . 
This experiment may be considered a landmark in the history 
of the Casimir effect, since it provided the first reliable exper- 
imental confirmation of this effect. One year later, using an 
atomic force microscope , Mohideen and Roy lfl6ll measured 
the Casimir force between a plate and a sphere with a better 
accuracy and established an agreement between experimen- 
tal data and theoretical predictions of less than a few percents 
(depending on the range of distances considered). The two 
precise experiments mentioned above have been followed by 
many others and an incomplete list of the modern series of ex- 
periments about the Casimir effect can be found in lfl7ll - ll26ll . 
For a detailed analysis comparing theory and experiments see 

We finish this subsection emphasizing that Casimir's orig- 
inal predictions were made for an extremely idealized situ- 
ation, namely: two perfectly conducting (flat) plates at zero 
temperature. Since the experimental accuracy achieved nowa- 
days is very high, any attempt to compare theory and exper- 
imental data must take into account more realistic boundary 
conditions. The most relevant ones are those that consider the 
finite conductivity of real metals and roughness of the surfaces 
involved. These conditions become more important as the dis- 
tance between the two bodies becomes smaller. Thermal ef- 
fects must also be considered. However, in principle, these 
effects become dominant compared with the vacuum contri- 
bution for large distances, where the forces are already very 
small. A great number of papers have been written on these 
topics since the analysis of most recent experiments require 
the consideration of real boundary conditions. For finite con- 
ductivity effects see Ref . ||30tl ; the simultaneous consideration 
of roughness and finite conductivity in the proximity for ap- 
proximation can be found in Ref. II31I1 and beyond PFA in Ref. 
II32I1 (see also references cited in the above ones). Concerning 
the present status of controversies about the thermal Casimir 
force see Ref. ^ 



B. The Casimir's approach 



The novelty of Casimir's original paper was not the predic- 
tion of an attractive force between neutral objects, once Lon- 
don had already explained the existence of a force between 



neutral but polarizable atoms, but the method employed by 
Casimir, which was based on the zero-point energy ot the elec- 
tromagnetic field. Proceeding with the canonical quantization 
of the electromagnetic field without sources in the Coulomb 
gauge we write the hamiltonian operator for the free radiation 
field as 



a=l k 



(3) 



where and ata are the creation and annihilation operators 
of a photon with momentum k and polarization a. The en- 
ergy of the field when it is in the vacuum state, or simply the 
vacuum energy, is then given by 



■Eo: 



^ 1 

EE 9^™k, 
k a=l ^ 



(4) 



which is also referred to as zero-point energy of the electro- 
magnetic field in free space. Hence, we see that even if we do 
not have any real photon in a given mode, this mode will still 
contribute to the energy of the field with j ftcOka and total vac- 
uum energy is then a divergent quantity given by an infinite 
sum over all possible modes. 

The presence of two parallel and perfectly conducting 
plates imposes on the electromagnetic field the following 
boundary conditions: 



Exn|,,,„,„. = 



(5) 



which modify the possible frequencies of the field modes. The 
Casimir energy is, then, defined as the difference between the 
vacuum energy with and without the material plates. How- 
ever, since in both situations the vacuum energy is a divergent 
quantity, we need to adopt a regularization prescription to give 
a physical meaning to such a difference. Therefore, a precise 
definition for the Casimir energy is given by 



■Ecas := lim 



E5^«k - El^^k 



, ka 



, ka 



//J 



(6) 



where subscript / means a regularized sum and that the fre- 
quencies are computed with the boundary conditions taken 
into account, subscript // means a regularized sum but with 
no boundary conditions at all and s stands for the regularizing 
parameter. This definition is well suited for plane geometries 
like that analyzed by Casimir in his original work. In more 
complex situations, like those involving spherical shells, there 
are some subtleties that are beyond the purposes of this intro- 
ductory article (the self-energy of a spherical shell depends on 
its radius while the self-energy of a pair of plates is indepen- 
dent of the distance between them). 

Observe that, in the previous definition, we eliminate the 
regularization prescription only after the subtraction is made. 
Of course, there are many different regularization methods. A 
quite simple but very efficient one is achieved by introducing 
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a high frequency cut off in the zero-point energy expression, 
as we shall see explicitly in the next section. This procedure 
can be physically justified if we note that the metallic plates 
become transparent in the high frequency limit so that the high 
frequency contributions are canceled out from equation (|6]l. 

Though the calculation of the Casimir pressure for the case 
of two parallel plates is very simple, its determination may 
become very involved for other geometries, as is the case, 
for instance, of a perfectly conducting spherical shell. Af- 
ter a couple of years of hard work and a "nightmare in Bessel 
functions", Boyer Issll computed for the first time the Casimir 
pressure inside a spherical shell. Surprisingly, he found a re- 
pulsive pressure, contrary to what Casimir had conjectured 
five years before when he proposed a very peculiar model for 
the stability of the electron ll34ll . Since then, Boyer's result 
has been confirmed and improved numerically by many au- 
thors, as for instance, by Davies in 1972 Issll . by Balian and 
Duplantier in 1978 ^ and also Milton in 1978 IstI], just to 
mention some old results. 

The Casimir effect is not a peculiarity of the electromag- 
netic field. It can be shown that any relativistic field under 
boundary conditions caused by material bodies or by a com- 
pactification of space dimensions has its zero-point energy 
modified. Nowadays, we denominate by Casimir effect any 
change in the vacuum energy of a quantum field due to any ex- 
ternal agent, from classical backgrounds and non-trivial topol- 
ogy to external fields or neighboring bodies. Detailed reviews 
of the Casimir effect can be found in iHMUllidliltl. 



C. A local approach 



In this section we present an alternative way of comput- 
ing the Casimir energy density or directly the Casimir pres- 
sure which makes use of a local quantity, namely, the energy- 
momentum tensor Recall that in classical electromagnetism 
the total force on a distribution of charges and currents can 
be computed integrating the Maxwell stress tensor through 
an appropriate closed surface containing the distribution. For 
simplicity, let us illustrate the method in a scalar field. The 
lagrangian density for a free scalar field is given by 



(7) 



The field equation and the corresponding Green function are 
given, respectively, by 

{-d^ + m^)isf{x) =0; (8) 



which, after symmetrization, can be written in the form 



(11) 



For our purposes, it is convenient to write the vacuum expec- 
tation value (VEV) of the energy-momentum tensor in terms 
of the above Green function as 



(0|7;,v W|0) = lim {d'dv+K^,) -8,v{d'ad°-+m^) G{x,x') . 

(12) 

In this context, the Casimir energy density is defined as 



Pcix) = (0|7oo(x)|0)bc - (0|rooW|0) 



Free i 



(13) 



where the subscript BC means that the VEV must be computed 
assuming that the field satisfies the appropriate boundary con- 
dition. Analogously, considering for instance the case of two 
parallel plates perpendicular to the O Z axis (the generaliza- 
tion for other configurations is straightforward) the Casimir 
force per unit area on one plate is given by 



jc = (o|7:t|o)-(o|7:T|o), 



(14) 



where superscripts + and — mean that we must evaluate (r,,) 
on both sides of the plate. In other words, the desired Casimir 
pressure on the plate is given by the discontinuity of (r,,) at 
the plate. Using equation (fT2b . {T„) can be computed by 

(0|r,,|0) = --lim(^^3^-^jG(.,.'). (15) 

Local methods are richer than global ones, since they pro- 
vide much more information about the system. Depending 
on the problem we are interested in, they are indeed neces- 
sary, as for instance in the study of radiative properties of an 
atom inside a cavity. However, with the purpose of comput- 
ing Casimir energies in simple situations, one may choose, 
for convenience, global methods. Previously, we presented 
only the global method introduced by Casimir, based on the 
zero-point energy of the quantized field, but there are many 
others, namely, the generalized zeta function method ll42ll and 
Schwinger's method ll43ll44ll . to mention just a few. 



II. EXPLICIT COMPUTATION OF THE CASIMIR FORCE 



(d'-m^)G{x,x') = -5{x-x'), 



(9) 



where, as usual, G(x,x') = i{0\T (^(^(x)(^{x')^ |0). 

Since the above lagrangian density does not depend explic- 
itly on X, Noether's Theorem leads naturally to the following 
energy-momentum tensor {d/jT^"^ = 0) 

dL 



(10) 



In this section, we show explicitly two ways of comput- 
ing the Casimir force per unit area in simple situations where 
plane surfaces are involved. We start with the global approach 
introduced by Casimir which is based on the zero-point en- 
ergy. Then, we give a second example where we use a local 
approach, based on the energy-momentum tensor. We finish 
this section by sketching some results concerning the Casimir 
effect for massive fields. 
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A. The electromagnetic Casimir effect between two parallel 
plates 



As our first example, let us consider the standard (QED) 
Casimir effect where the quantized electromagnetic field is 
constrained by two perfectly parallel conducting plates sepa- 
rated by a distance a. For convenience, let us suppose that one 
plate is located at z = 0, while the other is located at z = a. 
The quantum electromagnetic potential between the metallic 
plates in the Coulomb gau ge ( V • A = 0) which satisfies the 
appropriate BC is given by ||45I1 



A(p,z,f) = 



KckaL^ J 



1/2 



|a'*'(K,«)(K X z)sin 



3'2)(k,«) 



.n7t . fn%z\ ,K /nnz\ 
!— Ksm — z— cos 

ka \ a J k \ a / . 

+ h.c. , (16) 



where co(k,«) — ck ^ c [K^ + n^Tt^/a^] , with n being a 
non-negative integer and the prime in E' means that for « = 
an extra 1 /2 factor must be included in the normalization of 
the field modes. The non-regularized Casimir energy then 
reads 



he 

y 



n^Tl^ 



111 



he 



' adk, ^ 

yT' 



(27r)2y_„ 

Making the variable transformation + (nn/a)^ =: A, and in- 
troducing exponential cutoffs we get a regularized expression 
(in 1948 Casimir used a generic cutoff function). 



2jt 



1 

2 Jo 



IT 



Vdl 



dn 



r^^X-dX 



2k 




(17) 



Using that 



8e2 



'A' 

n=i 



-enn/a 



ae2 



1 



1 



e e^""/" - 1 



as well as the definition of Bernoulli's numbers. 



1 



- 1 



Lb,/—: 



we obtain 



E(fl) 



L2 
271 

271 




6(50-1)^^ + (l+2Z^0^ + ^(^ 



7=0 ) 

1 B4 



+ E^(-)'"'(«-2)(«-3)e^ 



e 

,«-4 



Using the well known values Bq=\, B\ 
and taking e ^ 0+, we obtain 



1 



and B4 



J_ 

30' 



L2 



= — ftc- 



7t"- 



1 



24x30 fl3 



As a consequence, the force per unit area acting on the plate 
at z = fl is given by 



F{a) 

L2 



1 d'Ec{c 
1? da 



TV- he 

' 240fl4 



-0,013 



1 



dyn 



[a/^imf cm^ 



where in the last step we substituted the numerical values 
of h and c in order to give an idea of the strength of the 
Casimir pressure. Observe that the Casimir force between 
the (conducting) plates is always attractive. For plates with 
lcm2 of area separated by l^m the modulus of this attrac- 
tive force is 0,013dyn. For this same separation, we have 
Pcas ~ ^O^^Patm, where Patm is the atmospheric pressure at 
sea level. Hence, for the idealized situation of two perfectly 
conducting plates and assuming = I err?, the modulus of 
the Casimir force would be « lO^^A^ for typical separations 
used in experiments. However, due to the finite conductivity 
of real metals, the Casimir forces measured in experiments are 
smaller than these values. 



B. The Casimir effect for a scalar field with Robin BC 



In order to illustrate the local method based on the energy- 
momentum tensor, we shall discuss the Casimir effect of a 
massless scalar field submitted to Robin BC at two parallel 
plates, which are defined as 



\bound. — P \bound. 

an 



(18) 



where, by assumption, P is a non-negative parameter How- 
ever, before computing the desired Casimir pressure, a few 
comments about Robin BC are in order 

First, we note that Robin BC interpolate continuously 
Dirichlet and Neumann ones. For p ^ we reobtain Dirichlet 
BC while for p ^ 00 we reobtain Neumann BC. Robin BC al- 
ready appear in classical electromagnetism, classical mechan- 
ics, wave, heat and Schrodinger equations [46] and even in the 
study of interpolating partition functions [47]. A nice realiza- 
tion of these conditions in the context of classical mechanics 
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can be obtained if we study a vibrating string with its extremes 
attached to elastic supports 1481 ? ]. Robin BC can also be used 
in a phenomenological model for a penetrable surfaces ll49ll . 
In fact, in the context of the plasma model, it can be shown 
that for frequencies much smaller than the plasma frequency 
(co ^ ©/>) the parameter P plays the role of the plasma wave- 
length [? ]. In the context of OFT this kind of BC appeared 
more than two decades ago ifsoHsill . Recently, they have been 
discussed in a variety of contexts, as in the AdS/CFT corre- 
spondence ll52ll . in the discussion of upper bounds for the ratio 
entrop y/en ergy in confined systems ||53I1 . in the static Casimir 
effect fsi!], in the heat kernel expansion llssi Is^ IstIi a nd in 
the one-loop renormalization of the A,(|)^ theory llssl 15911 . As 
we will see, Robin BC can give rise to restoring forces in the 
static Casimir effect ifsill . 

Consider a massless scalar field submitted to Robin BC on 
two parallel plates: 



we have 



az oz 



(19) 



where we assume Pi (P2) > 0. For convenience, we write 

d^ka 



G{x,x) = J 



(27r)3 



„'A:a(A-x')' 



g(z,z';k^,co), (20) 



with a = 0, 1,2, = diag{— I, and the reduced 
Green function §(z,z';kx,CO) satisfies 



Kz,z';kx,co) = -5{z-z'), 



(21) 



where A,^ = co^ — andg(z,z';kx,Co) is submitted to the fol- 
lowing boundary conditions (for simplicity, we shall not write 
kx, CO in the argument of g): 

d 



g{Q,z')^^i^^g{Q,z') 



?{a,z') 



-h^g{a,z') 
az 



It is not difficult to see that g{z,z') can be written as 

A(z')(sinA,z + piA,cosA,z), z<z' 



^iz,z')-- 



B{z') (sinX(z-fl)-P2XcosA.(z-fl)), z>z' 



For points inside the plates, < < a, the final expression 
for the reduced Green function is given by 



,RR 



(.z,z') 



J.^'>"(<":>-") . 
T2 



72 



(22) 



where y,- = ^ -jj'^ (' = 1 , 2). Outside the plates, with z,z' > a, 
we have: 



iz,z) 



2iX 

Defining f''^ such that 



J2 



{2n) 2k ^ 



,-iA,(z<-a) 



(23) 



(24) 



The Casimir force per unit area on the plate at x = a is given 
by the discontinuity in (f^^): 



T = 



d^k I dm 



{2k) 2% 

After a straightforward calculation, it can be shown that 



/ l+Pl^/2« \/ l+p24/2« \ 



(26) 

Depending on the values of parameters Pi and P2, restoring 
Casimir forces may arise, as shown in Figure 2 by the dotted 
line and the thin solid line. 




FIG. 2: Casimir pressure, conveniently multiplied by a , as a func- 
tion of a for various values of parameters (3i and (32. 

The particular cases of Dirichlet-Dirichlet, Neumann- 
Neumann and Dirichlet-Neumann BC can be reobtained if 
we take, respectively, pi = p2 = 0, Pi = p2 °o and pi = 
; P2 ^ °o. For the first two cases, we obtain 

^^^(a) = f^^ia) = --i^ ^ . (27) 

i2K^a^ Jo e'^—X 

Using the integral representation Jq d^ — - — C,R{s)r{s), 

where II^r is the Riemann zeta function, we get half the elec- 
tromagnetic result, namely. 



For the case of mixed BC, we get 



1 



TV- he 1 
' 480 ^ 



d^ 



(28) 



(29) 
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Using in the previous equation the integral representation 

T+T ^ ^ 2'^'*)r(i)^^(i) we get a repulsive pres- 
sure (equal to half of Boyer's result |[6(1| obtained for the elec- 
tromagnetic field constrained by a perfectly conducting plate 
parallel to an infinitely permeable one), 



7 %^hc 1 



(30) 
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C. The Casimir effect for massive particles 



In this subsection, we sketch briefly some results concern- 
ing massive fields, just to get some feeling about what kind 
of influence the mass of a field may have in the Casimir ef- 
fect. Firstly, let us consider a massive scalar field submit- 
ted to Dirichlet BC in two parallel plates, as before. In this 
case, the allowed frequencies for the field modes are given 



by cOk = c 



1/2 



which lead, after we use the 



Casimir method explained previously, to the following result 
for the Casimir energy per unit area (see, for instance, Ref. 
B) 



1 , , 



L-^ K2[2amn) , (31) 



where is a modified Bessel function, m is the mass of the 
field and a is the distance between the plates, as usual. The 
limit of small mass, am ^ 1, is easily obtained and yields 



what kind of BC can be chosen for this field. The point is 
that Dirac equation is a first order equation, so that if we want 
non-trivial solutions, we can not impose that the field satis- 
fies Dirichlet BC at two parallel plates, for instance. The most 
appropriate BC for fermions is borrowed from the so called 
MIT bag model for hadrons ll62ll . which basically states that 
there is no flux of fermions through the boundary (the normal 
component of the fermionic current must vanish at the bound- 
ary). The Casimir energy per unit area for a massive fermionic 
field submitted to MIT BC at two parallel plates was first com- 
puted by Mamayev and Trunov |63| (the massless fermionic 
field was first computed by Johnson in 1975 ll64ll ) 



1 



1 



n^a^ 



J ma 



1 



^ — t7ja 



2^ 



The small and large mass limits are given, respectively, by 



^2;/ («,'«) 



771^ 



2880fl3 24a2 
3{mayf_ 



^ — 2ma 



(34) 



Since the first correction to the zero mass result has an oppo- 
site sign, also for a fermionic field small masses diminish the 
Casimir effect. In the large mass limit, ma oo, we have a 
behaviour analogous to that of the scalar field, namely, an ex- 
ponential decay with ma (again this happens due to the plane 
geometry). We finish this section with an important observa- 
tion: even a particle as light as the electron has a completely 
negligible Casimir effect. 



1 



■Ec{a,m) ; 



71" 



1440a^ 



m 

96a 



(32) 



As expected, the zero mass limit coincides with our previous 
result ( |28] ) (after the force per unit area is computed). Observ- 
ing the sign of the first correction on the right hand side of last 
equation we conclude that for small masses the Casimir effect 
is weakened. 

On the other hand, in the hmit of large mass, am S> 1, it can 
be shown that 



'£c{a,m) 



m^ / K 
l6K^a \ma 



1/2 



(33) 



Note that the Casimir effect disappears for m °o, since in 
this limit there are no quantum fluctuations for the field any- 
more. An exponential decay with ma is related to the plane 
geometry. Other geometries may give rise to power law de- 
cays when ma °°. However, in some cases, the behaviour 
of the Casimir force with ma may be quite unexpected. The 
Casimir force may increase with ma before it decreases mono- 
tonically to zero as ma °o (this happens, for instance, when 
Robin BC are imposed on a massive scalar field at two parallel 
plates ieill ). 

In the case of massive fermionic fields, an analogous be- 
haviour is found. However, some care must be taken when 
computing the Casimir energy density for fermions, regarding 



III. MISCELLANY 



In this section we shall briefly present a couple of topics 
which are in some way connected to the Casimir effect and 
that have been considered by our research group in the last 
years. For obvious reasons, we will not be able to touch all 
the topics we have been interested in, so that we had to choose 
only a few of them. We first discuss how the Casimir effect 
of a charged field can be influenced by an external magnetic 
field. Then, we show how the constitutive equations asso- 
ciated to the Dirac quantum vacuum can be affected by the 
presence of material plates. Finally, we consider the so called 
dynamical Casimir effect. 



A. Casimir effect under an external magnetic field 



The Casimir effect which is observed experimentally is that 
associated to the photon field, which is a massless field. As 
we mentioned previously, even the electron field already ex- 
hibits a completely unmeasurable effect. With the purpose 
(and hope) of enhancing the Casimir effect of electrons and 
positrons we considered the influence of an external elec- 
tromagnetic field on their Casimir effect. Since in this case 
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we have charged fields, we wondered if the virtual electron- 
positron pairs which are continuously created and destroyed 
from the Dirac vacuum would respond in such a way that 
the corresponding Casimir energy would be greatly amplified. 
This problem was considered for the first time in 1998 ll65ll 
(see also Ref. lH^). The influence of an external magnetic 
field on the Casimir effect of a charged scalar field was con- 
sidered in Ref. |f67j| and recently, the influence of a magnetic 
field on the fermionic Casimir effect was considered with the 
more appropriate MIT BC ll68ll . 

For simplicity, let us consider a massive fermion field un- 
der anti-periodic BC (in the OZ direction) in the presence 
of a constant and uniform magnetic field in this same direc- 
tion. After a lengthy but straightforward calculation, it can be 
shown that the Casimir energy per unit area is given by ll66ll 

E(fl,B) 2(flm)2 ~ (-1)"-' , 



^2 



4n\f^^^ 'Jo \ a 



where we introduced the Langevin function: L(^) = coth^ — 
1 In the strong field limit, we have 



E(fl,B) 



eBm ~ (-1)"-' 
^ > Ki(amn) 



(35) 



In this limit we can still analyze two distinct situations, 
namely, the small mass limit {ma <C 1) and the large mass 
limit {ma 3> 1). Considering distances between the plates typ- 
ical of Casimir experiments (a « 1 iim), we have in the former 



9c{a,B) 
Pc(a,0) 



10" 



Tesla 



ma <IC 1 , 



(36) 



while in the latter case. 



9c{a,B) 
Pc(fl,0) 



lO-'Ox 



B 



Tesla 



ma 3> 1 . 



(37) 



In the above equations Pt(fl,Z?) and Pc(a,0) are the Casimir 
energy density under the influence of the magnetic field and 
without it, respectively. Observe that for the Casimir ef- 
fect of electrons and positrons, which must be treated as the 
large mass limit described previously, huge magnetic fields 
are needed in order to enhance the effect (far beyond accessi- 
ble fields in the laboratory). In other words, we have shown 
that though the Casimir effect of a charged fermionic field can 
indeed be altered by an external magnetic field, the universal 
constants conspired in such a way that this influence turns out 
to be negligible and without any chance of a direct measure- 
ment at the laboratory. 



B. Magnetic permeability of the constrained Dirac vacuum 



In contrast to the classical vacuum, the quantum vacuum 
is far from being an empty space, inert and insensible to any 
external influence. It behaves like a macroscopic medium, in 
the sense that it responds to external agents, as for example 
electromagnetic fields or the presence of material plates. As 
previously discussed, recall that the (standard) Casimir effect 
is nothing but the energy shift of the vacuum state of the field 
caused by the presence of parallel plates. There are many 
other fascinating phenomena associated to the quantum vac- 
uum, namely, the particle creation produced by the application 
of an electric field [[69tl, the birefringence of the QED vacuum 
under an external magnetic field iFtoI ItTIi and the Scharnhorst 
effect ll72i ItsIi . to mention just a few. This last effect pre- 
dicts that the velocity of light propagating perpendicularly to 
two perfectly conducting parallel plates which impose (by as- 
sumption) BC only on the radiation field is slightly altered by 
the presence of the plates. The expected relative variation in 
the velocity of light for typical values of possible experiments 
is so tiny (Ac/c w 10^^^) that this effect has not been con- 
firmed yet. Depending on the nature of the material plates the 
velocity of light propagating perpendicular to the plates is ex- 
pected to diminish [74]. The Scharnhorst effect inside a cavity 
was considered in [75j]. 

The negligible change in the velocity of light predicted by 
Scharnhorst may be connected with the fact that the effect that 
bears his name is a two-loop QED effect, since the classical 
field of a traveling light wave interacts with the radiation field 
only through the fermionic loop. 

With the purpose of estimating a change in the constitutive 
equations of the quantum vacuum at the one-loop level, we 
were led to consider the fermionic field submitted to some BC. 
The magnetic permeability /j of the constrained Dirac vacuum 
was computed by the first time in Ref. Jt^, but with the non- 
realistic anti-periodic BC. In this case, the result found for the 
relative change in the permeability, A^u := — 1, was also neg- 
ligible (an analogous calculation has also been made in the 
context of scalar QED IttII ). However, when the more real- 
istic MIT BC are imposed on the Dirac field at two parallel 
plates things change drastically. In this case, it can be shown 
that the magnetic permeability of the constrained Dirac vac- 
uum is given by IITSIl 



2 



/u{ma) 



12jt2 



where 



Hima) 



1 '^''(x2- 1)3/2 



ln<^ 1 



-H{ma) , 



For confining distances of the order of 0, 1/jot, we have 



A^l■.= ^l-\^\Q 



-9 



(38) 



(39) 



The previous value is comparable to the magnetic permeabil- 
ity of Hydrogen and Nitrogen at room temperature and atmo- 
spheric pressure. Hence, an experimental verification of this 
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result seems to be not unfeasible (we will come back to this 
point in the final remarks). 



C. The Dynamical Casimir effect 



As our last topic, we shall briefly discuss the so called dy- 
namical Casimir effect, which consists, as the name suggests, 
of the consideration of a quantum field in the presence of 
moving boundaries. Basically, the coupling between vacuum 
fluctuations and a moving boundary may give rise to dissi- 
pative forces acting on the boundary as well as to a particle 
creation phenomenon. In some sense, these phenomena were 
expected. Recall that the static Casimir force is a fluctuat- 
ing quantity ll79ll and hence, using general arguments related 
to the fluctuation-dissipation theorem Isoll dissipative forces 
on moving boundaries are expected. Further, using arguments 
of energy conservation we are led to creation of real parti- 
cles (photons, if we are considering the electromagnetic field 
Isill ). For the above reasons, this topic is sometimes referred 
to as radiation reaction force on moving boundaries. 

After Schwinger's suggestion that the phenomenon of 
sonoluminescence could be explained by the dynamical 
Casimir effect [fsj] (name coined by himself), a lot of work 
has been done on this subject. However, it was shown a few 
years later that this was not the case (see and references 
therein for more details). 

The dynamical Casimir effect already shows up in the case 
of one (moving) mirror ||83[ [84l [ssll . However, oscillating 
cavities whose walls perform vibrations in parametric reso- 
nance with a certain unperturbed field eigenfrequency may 
greatly enhance the effect IH, [13, HI, H^- Recendy, a one 
dimensional oscillating cavity with walls of different nature 
was considered The dynamical Casimir effect has also 
been analyzed for a variety of three-dimensional geometries, 
including parallel plane plates ll9lll . cylindrical waveguides 
ll92ll . and rectangular jpst], cylindrical ll94ll and spherical cavi- 
ties Igsll . For a review concerning classical and quantum phe- 
nomena in cavities with moving boundaries see Dodonov ||96|] 
and for a variety of topics on non-stationary Casimir effect 
including perspectives of its experimental verification see the 
special issue ll97ll . 

In this section, we shall discuss the force exerted by the 
quantum fluctuations of a massless scalar field on one moving 
boundary as well as the particle creation phenomenon in a un- 
usual example in 1h-1 dimensions, where the field satisfies a 
Robin BC at the moving boundary IHIll. We shafl follow 
throughout this paper the perturbative method introduced by 
Ford and Vilenkin lf84l] . This method was also applied suc- 
cessfully to the case of the electromagnetic field under the 
influence of one moving (perfectly) conducting plate llssll as 
well as an oscillating cavity formed by two parallel (perfectly) 
conducting plates 19111 . 

Let us then consider a massless scalar field (|) in 1 + 1 in 
the presence of one moving boundary which imposes on the 
field a Robin BC at one moving boundary when observed from 
a co-moving inertial frame. By assumption, the movement 



of the boundary is prescribed, non-relativistic and of small 
amplitude (5^(f) is the position of the plate at a generic in- 
stant f). Last assumptions may be stated mathematically by 
|54'(0I << c and |8^(f)| << c/coo , where coq corresponds 
to the (main) mechanical frequency. Therefore, we must solve 
the following equation: d^(^{x,t) = 0, with the field satisfying 
a Robin BC at the moving boundary given by 



*l'(^,0l-c=5?(0 = p'l'(-^,0lx=8<7(r) • (40) 



where p is a non-negative parameter and the previous condi- 
tion was already written in the laboratory frame. We are ne- 
glecting terms of the order {dq^/c^). The particular cases of 
Dirichlet and Neumann BC are reobtained by making p = 
and P ^ oo, respectively. The dissipative forces for these par- 
ticular cases were studied at zero temperature as well as at 
fini te te mperature and also with the field in a coherent state 
in jlOOl l (dissipative forces on a perfectly conducting mov- 
ing plate caused by the vacuum fluctuations of the electro- 
magnetic field at zero and non-zero temperature were studied 
in lIlOIll ). The perturbative approach introduced by Ford and 
Vilenkin llsill consists in writing 



(j)(x,f) =(|)o(jc,r) + 5(|)(x,f) 



(41) 



where (|)o(x,f) is the field submitted to a Robin BC at a static 
boundary fixed at the origin and 5(|)(.Y,r) is the first order con- 
tribution due to the movement of the boundary. The total force 
on the moving boundary may be computed with the aid of the 
corresponding energy-momentum tensor, namely. 



t,sq-mo) 



(42) 



where superscripts + and — mean that we must compute the 
energy-momentum tensor on both sides of the moving bound- 
ary. It is convenient to work with time Fourier transforms. 
The susceptibility x(co) is defined in the Fourier space by 



(43) 



where 5T (co) and dQ{(i)) are the Fourier transformations of 
dF{t) and 5^(f), respectively. It is illuminating to compute 
the total work done by the vacuum fluctuations on the moving 
boundary. It is straightforward to show that 



8F{t)8q{t)dt 



K Jo 



t/coco /mx(a)) |52((X)) p . (44) 



Note that only the imaginary part of x{g)) appears in the pre- 
vious equation. It is responsible for the dissipative effects and 
hence it is closely related with the total energy converted into 
real particles. On the other hand, the real part of when it 
exists, does not contribute to the total work and hence it is not 
related to particle creation, but to dispersive effects. For the 
particular cases of Dirichlet or Neumann BC it can be shown 
that the susceptibility is purely imaginary and given by 



X^(co)=x''(co) = / 



6kc^ 



(45) 



10 



etal. 



which imphes 



bF{t) 



bq{t) 



(46) 



Since /ot%(co) > 0, for these cases the vacuum fluctuations 
are always dissipating energy from the moving boundary. 

However, for Robin BC an interesting thing happens. It 
can be shown that acquires also a real part, which gives 
rise to a dispersive force acting on the moving boundary. The 
explicit expressions of 1{_ex{(£)) and lmx{(£>) can be found 
in iosl] . but the general behaviour of them as functions of 
CO is shown in Figure [5] For convenience, we normalize 
these quantities dividing them by the value irnXoi^), where 
the subscript D means that /»zx((o) must be computed with 
Dirichlet BC. 





/ Ini x(<a) / Im Xn(co) 

J / 




1 1 Re / Im x„(ro) 




i 



FIG. 3: Imaginary and real parts of x(co) with Robin BC appropri- 
ately normalized by the value of Im%((s3) for the Dirichlet BC. 

Now, let us discuss briefly the particle creation phe- 
nomenon under Robin BC. Here, we shall consider a semi- 
infinite slab extending from — oo to ?>q{t) following as before a 
prescribed non-relativistic motion which imposes on the field 
Robin BC at 5^(f). It can be shown that the corresponding 
spectral distribution is given by ll99ll 



dN 
d(£) 



(CO) 



4co 



-dti)' [8e(co-co')] 



■CO 



l-p2(^co' 



(47) 

As an explicit example, let us consider the particular motion 

hq{t) = S^oe"'''''^ cos(cOof ) , (48) 

where, by assumption, cooT 1 (this is made in order to sin- 
gle out the effect of a given Fourier component of the motion). 
For this case, we obtain the following spectral distribution 



dN 
d(£) 



(co)= (8^o)^7'tt)(tt)o-tt>)x 

[1-pMcoo-co)]^ 

''(l+p2co2)(l + p2(cOo-CO)2) 



©(coo-co). (49) 



The spectral distributions for the particular cases of Dirichlet 
or Neumann BC can be easily reobtained by making simply 



P = an d P oa, respectively. The results coincide and are 
given by ll02ll (we are making c = 1 ) 



d(0 



(CO) 



dN^ 
d(£) 



(co) = (8^o)^7'tt>(coo - co)0(coo - co) 



(50) 

A simple inspection in ( |49l ) shows that, due to the presence of 
the Heaviside step function, only the field modes with eigen- 
frequencies smaller than the mechanical frequency are ex- 
cited. Further, the number of particles created per unit fre- 
quency when Robin BC are used is always smaller than the 
number of particles created per unit frequency when Dirichlet 
(or Neumann) BC are employed. 




FIG. 4: Spectral distributions of created particles for: Dirichlet and 
Neumann BC (dashed line) and for some interpolating values of the 
parameter (3 (dotted and solid lines). 

Figure |4] shows the spectral distribution for different val- 
ues of the parameter p, including P = (dashed line), which 
corresponds to Dirichlet or Neumann BC. Note that for inter- 
polating values of P particle creation is always smaller than 
for p = (dotted line). Depending on the value of p, particle 
creation can be largely suppressed (solid line). 



IV. FINAL REMARKS 



In the last decades there has been a substantial increase in 
the study of the Casimir effect and related topics. It is remark- 
able that this fascinating effect, considered nowadays as a fun- 
damental one in QFT, was born in connection with colloidal 
chemistry, an essentially experimental science. As we men- 
tioned previously, the novelty of Casimir's seminal work fl]] 
(see also |j^) was the technique employed by him to compute 
forces bet ween neutral bodies as is emphasized by Itzykson 
and Zuber ifiol: 



By considering various types of bodies influenc- 
ing the vacuum configuration we may give an 
interesting interpretation of the forces acting on 
them. 
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However, the attractive or repulsive character of the Casimir 
force can not be anticipated. It depends on the specific bound- 
ary conditions, the number of space-time dimensions, the na- 
ture of the field (bosonic or fermionic), etc. The "mystery" 
of the Casimir effect has intrigued even proeminent physicists 
such as Julian Schwinger, as can be seen in his own statement 
0: 

... one of the least intuitive consequences of quan- 
tum electrodynamics. 

There is no doubt nowadays about the existence of the 
(static) Casimir effect, thanks to the vast list of accurate ex- 
periments that have been made during the last ten years. It is 
worth emphasizing that a rigorous comparison between theory 
and experimental data can be achieved only if the effects of 
temperature and more realistic BC are considered. In princi- 
ple, the former are important in comparison with the vacuum 
contribution for large distances, while the latter can not be 
neglected for short distances. Typical ranges investigated in 
Casimir experiments are form 0. 1/jm to 1 .0/jm and, to have an 
idea of a typical plasma wavelengths (the plasma wavelength 
is closely related to the penetration depth), recall that for Au 
we have Xp k, \3>6nm. 

The Casimir effect has become an extremely active area 
of research from both theoretical and experimental points of 
view and its importance lies far beyond the context of QED. 
This is due to its interdisciplinary character, which makes this 
effect find applications in quantum field theory (bag model, 
for instance), cavity QED, atomic and molecular physics, 
mathematical methods in QFT (development of new regular- 
ization and renormalization schemes), fixing new constraints 
in hypothetical forces, nanotechnology (nanomachines oper- 
ated by Casimir forces), condensed matter physics, gravitation 
and cosmology, models with compactified extra-dimensions, 
etc. 

In this work, we considered quantum fields interacting only 
with classical boundaries. Besides, these interactions were 
described by highly idealized BC. Apart from this kind of in- 
teraction, there was no other interaction present. However, 
the fields in nature are interacting fields, like those in QED, 
etc. Hence, we could ask what are the first corrections to the 
Casimir effect when we consider interacting fields. In princi- 
ple, they are extremely small. In fact, for the case of QED, the 



first radiative correction to the Casimir energy density (con- 
sidering that the conducting plates impose BC onl y on the ra- 
diation field) was firstly computed by Bordag et al lll04ll and is 

given by £^'^(0,0) — T.'^'' ^^^^2~cf'' where Xc is the Comp- 

ton wavelength of the electron and E^"' (a) is the zeroth order 
contribution to the Casimir energy density. As we see, at least 
for QED, radiative corrections to the Casimir effect are exper- 
imentally irrelevant. However, they might be relevant in the 
bag model, where for quarks Xc ~ a and also the quark prop- 
agators must be considered submitted to the bag BC ||9]- Be- 
sides, the study of radiative corrections to the Casimir effect 
provide a good laboratory for testing the validity of idealized 
BC in higher order of perturbation theory. 

Concerning the dynamical Casimir effect, the big challenge 
is to conceive an experiment which will be able to detect real 
photons created by moving boundaries or by an equivalent 
physical system that simulates rapid motion of a boundary. An 
ingenious proposal of an ex periment has been made recently 
by the Padova's group llOSll . There are, of course, many other 
interesting aspects of the dynamical Casimir effect that has 
been studied, as quant um de coherence IIO6II . mass correction 
of the moving mirrors llOTll . etc. (see also the reviews llOSll ). 

As a final comment, we would like to mention that surpris- 
ing results have been obtained when a deformed quantum field 
theory is considered in connection with the Casimir effect. It 
seems that the simultaneous assumptions of deformation and 
boundary conditions lead to a new mec hanis m of creation of 
real particles even in a static situation ll09ll . Of course, the 
Casimir energy density is also modified by the deformation 
il lOll . Quantum field theories with different space-time sym- 
metries, other than those governed by the usual Poinc are al - 
gebra (as for example the K-deformed Poincare algebra ill 1 ill ) 
may give rise to a modified dispersion relation, a desirable fea- 
ture in some tentative models for solving recent astrophysical 
paradoxes. 
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